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Abstract 

It is shown that nonlinear terms in equations of gravitons on the back- 
ground of curved space-time of the expanding Universe can solve the prob- 
lem of the negative square of the effective mass formally arising in linear 
approximation for gravitons. Similar to well known spontaneous breaking 
of symmetry in Goldstone model one must take another vacuum so that 
nonzero vacuum expectation value of the quantized graviton field leads to 
change of spectrum for gravitons. There appears two graviton fields, one 
with the positive mass, another with the zero mass. Energy density and 
the density of particles created by gravitation of the expanding Universe 
are calculated for some special cases of the scale factor. Numerical of 
result are obtained for the dust universe case. 

Introduction 

There is a problem in quantum theory of gravitons created from vacuum in the 
expanding Universe with nonzero scalar curvature R (inflation, dust, etc.) con- 
cerning the long wave graviton modes. In linearized theory of quantum gravitons 
in curved space-time of the isotropic homogeneous Universe one obtains after 
separation of variables in the wave equation the equation for the function de- 
pendent only on time. This equation can be understood as equation in flat 
stationary metric with time dependent mass. It occurs that for long waves this 
effective mass squared is negative. All this occurs due to conformal noninvari- 
ance of the graviton theory for nonzero R leading to tachyonic behaviour of long 
waves modes. 

In some papers pQ (see and references there) it was proposed to consider 
these modes as classical excitations of the field growing in time, so that one must 
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quantize only modes with momentum with the square larger than the negative 
square of the effective mass. However one knows from the quantum field theory 
that tachyonic behaviour disappears if one takes into account nonlinear terms 
neglected in linearized theory. This is typical in theories of spontaneous breaking 
of symmetry due to redefinition of the vacuum leading to its noninvariance to 
this or that transformation of the Lagrangian. In quantum theory based on a 
new vacuum one gets new masses for the redefined quantum field so that there 
is no negative mass square. 

In this paper the analogous program is made for gravitons. It occurs that 
if one is going from the linear theory of gravitons taking into account the next 
order of nonlinearity one gets the redefinition of vacuum solving the problem 
for long wave gravitons. In the result one gets gravitons with zero and positive 
effective mass. In the end of the paper the expressions for the particle density 
and the energy density are obtained for gravitons created in expanding Universe 
with metric which has some special dependence of the scale factor on time. 



1 Getting the graviton equation from Einstein 
equation 

Einstein equations in presence of matter have the form 

Rik - -j9ikR = nT lk (1) 

or 

k = <n - ±#r). (2) 

Let us consider the case when matter is homogeneous isotropic liquid filling 
the Universe. Then 

T ik = (e + p)uiU k - gikP, (3) 

where u, is the four velocity, p - the pressure and e - the energy density of the 
liquid. 

The problem of creation of gravitons in the early Universe was discussed 
in literature with gravitons considered as quantized small term in the metrical 
tensor. Due to absence of exact quantum gravity usually one deals with lin- 
earized analogy with quantization of other quantum fields. First let us obtain 
equations for classical small perturbations of the metrical tensor and then do 
quantization. Consider the graviton perturbations-the gravitational waves as 
small term added to the background metric (4) So there are 

(°) 

9ik = 9 ik +h lk (4) 

(°) 

if hik = the 9 ik is the solution of Einstein's equation of the form 

(°) ■ (°) ■ 1 • (°) 

Rl = <Tl--Sl T). (5) 
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Let us go from the up to low indices and vice verse by using the background 

(°) (°) . 

metric 9 ik'- h l k = g m h n k and expand equations J5| in a series in h\: 

(°) ■ (°) • 1 - (°) 1 ■ 

R I + SRI = k(T \- -S\ T +STI - 2 S i § n (6) 

from which due to ^ the perturbations hik satisfy equations 

SRI = k(STI - 1 -S18T). (7) 

Using the notation (1 + h)^ 1 for the matrix inverse to (1 + h) with small h\ 
(small in the sense that all eigenvalues of the matrix (1 + h) are smaller than 
the unit)one obtains 

(i+/i)- i i = 4-^i+w---- (8) 

Write the Ricci tensor and the curvature tensor as 

R{ = (1 + h)-%{-ti ( Rl + \{l + h)- ll v {h% + ft£jf - 4* - <f )+ 

+ h)- ll v (l + hrZihUhf'*' - (2hi, - h{ n )- 
■ (h?>* + hf - h^') - 2h l Uhf - hf ; "')))+ R I (9) 
R lk =Rik + h)- 1 ],^.^ + - hfi tl ~ hl Mi )+ 

+ h)- l \,{i + hyZKuht - ( 2h U - h U- 

■ Ki + K k - hf) - 2h l k \ n (h^ - (10) 
R = R\=R +(1 + h)- l \,{-hi R 2 + (1 + h)- l \,{h l -f - h]f)+ 

+ h)- il v {i + hrZishiX 11 ' - Hnhf- 

(°) 

Here ";" means the covariant derivative in background metric 9 ik- Consid- 
ering in (|10p only first degree in hik one obtains the linearized equations for hik 
in (TTJ) as 

h%, n + K, k - K Mn - hl^- { °g\k (hi- - h££) + h ik R= -2nS T lk ■ (12) 

Let us consider the background as homogeneous isotropic nonstastionary 
space-time 

ds 2 = dt 2 - a 2 (t)dl 2 , (13) 
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where 77 is the conformal time. Here the Latin indices take the values 0, 1, 2, 3 
and the Greek - 1,2,3. Then write for the scalar curvature and the Ricci tensor 

5R=(1 + h)- lk k ,((l + fc)" 1 !,^ - h? 1 '?) - ^<o, - ^< + 



-2hlr n h%' k ' - 4 n hf< n ') + ^(1 + hr^iShlXo - hlohQ). 



SR°o = "2^(1 + h)- ll lf (hi Q ,o + ~ ^(1 + ^Xo^o)- 

^ = ^K 1 + ^UCo ~ Co) + 4^(1 + jo -1 "*- 

•C4a^O-<o(<TI-^))- 

5B% = (1 + ^"(-(l + ft)- x !,((/»fS + hf\ - ti?\ - h$)+ 
ih#* - hp n ')) - ^(h$Xo - 2<o^,o))) - ^^'0,0 " ^,0 + f ^) 



-^(l + Z^X' (14) 

where e = ±1, for the closed, open and flat Universe. The sign is used for 
the covariant derivative in space part of the metric and " , 0" or comma for the 
derivative in conformal time rj. Metric g lk is defined up to arbitrary coordinate 
transformations so one can put the condition. 
Solutions ([6]) for hik can be written as 

hi = Si + V£ + B{, 

where S l k ,V£,B l k are irreducible scalar, vector and tensor components of the 
tensor satisfying the conditions [2]: 

(°) 1 (°) -v (°) i. 

g lk S tk = S^0, g lk V lk = 0, 9 ik B lk = 0. 

Considering only the gravitational waves exclude the scalar and vector parts 
by putting the gauge conditions 

h = 0, fcj,. = 0. (15) 

In this case the linearized equations for perturbations h\ take the form 

^,o,o + ^f,o + + h% = 0. (16) 
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Go from the variables h\ to new variables fi\ — a(r))h\ and make the con- 
formal transformation 

9tk = gik/a 2 (rj) 

Then one obtains the equation for the field with spin 2 in Minkowsky flat 
space in some external effective field 

M?,o,o + ^-^ = 0. (17) 
After separation of variables and Fourier representation jj/g (x) 

l$(x) = J ^k(jg % {r,)e^a% + g^e^ 'af) (18) 
one obtains for the time dependent gt(f]) function the equation 

gl(v) + (k 2 -^)g n (v) = 0, (19) 
which formally has the negative square of the effective mass m 2 a 2 = —a" /a. 

2 Third order equations for gravitational waves 

Let us consider the right hand side of Einstein's equations. For {3| 

STj. = (5e + dpjrfiik - 6* k 6p + (e + p)(Su l u k + u l 5u k + 5u l 5u k ). (20) 

(°) . . i 

One has for the four velocities Ui and u j the conditions gi k u l u = 1 and 

( n } (0) * (0) fc 1 a 

g ik u 1 u * = l. So 

2 ( u k 6u k + ( °g\ k 5u l 5u k = 0. (21) 

Note that 8ui due to constraints (fTSf can depend only on h 3 k h k H ... so 8u a 6u^ 
depends on the squares of these terms. But we shall neglect the fourth and 
higher orders. In synhronous reference system u° = l/a,u a = so from (|21[) 
one has 5u° — and 

8R% = ~K{5e + 36p), 5R^ = ^K6%{6p - Se), 5R% = an{e + p)8u a . (22) 

Consider the case of fiat space e = 0. Then h%„ = h% „, hPj 1 — —h^' 1 where 

Greek indices are put up and below by use of the Minkowsky metric. One can 
look on eqs. ([22]) as on Euler Lagrange equations for the fields hi k . Then due 
to constraints (|15|) up to terms of the divergence form one can obtain that not 
only hy = but h k h^ k = 0, ... so that (H]) can be transformed to 

5R$ = (1 + h)-^(hl + h)-%{~h^ - 
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+|(1 + h)- l \,{l + hrZih^hf'' - h<t ln rf n ' + 2h% n h( ln '))- 

-^(1 + ^)" 1 !'(^,o^o - 2ft|>jg l0 ))- (23) 
From (I23p follows that one can take instead of ([7]) the equations 

(l + h)Z6R% = ±K(l + h)}(5p-5 E ). (24) 
Multiply the last equation on "—2a 2 ", then from (|23[) . ([24)1 one obtains 



hlo, a + ^h% + ^(1 + + 



a a 



+ /O^^Xo - 2^o^,o) + (1 + fc) _1 r(^;r + O" 

+ h)- l \,{l + h)-Z(KX ,a ~ Knh a / + 2h inK y ) = 

= a 2 K(l + h)^(Se-Sp). (25) 
Consider first three orders in h\ in equations (I25p . 



2a' 1 



^,o,o + — h% + -(Si, h\, + h l k hf,)(h% hl - 2hf fi h l fit0 + 2h% + 2h\ ; ; 



or 



= [8f + h%)(a 2 K (Se - Sp) + ^(h\, - h l k 4)h( ), (26) 



h a 4- ^-h a 4- h a ' L h a h l h a h l ' n \h n ' a h l h l h l ' a 

"73,0,0 t n 0, ' n i3,l ~ n l,O n 0,O ~ a l,n n ~2 1 n <P n l 1,0 



-5%(a<K(Se-Sp) + -hZ,hl )+ 



+ l -h\,{{2hi^ a - 2h?'%, n - h{ n h a n + ^h\f - hlXfi + K,oKo+ 
+ — 5%h?hi fi - —h%)h{ ) = h%a 2 K (Se - Sp) = 0. (27) 



So 



6%a 2 K(5p-5e) = hfah^o + hf, n h l f + + ti'Kf + ^~>&<o (28) 
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Putting away the divergence of h^h 1 ^ + ^hj^h^ ^ + and taking into 

account for fixed nonzero components of the tensor h l k the condition h\,hfh l n = 
after simple transformations one obtains 



2a' 



+h\,{\hi n h«> n + h i :>«hi P - h^xr + h i n h^) = o. (29) 

Now let us go from variables h\ to variables /4 = a ( 7 ?)^ l l an d make the 
conformal transformation 

9ik = 9ik/a 2 (i]). 

Then we obtain the equation in flat Minkowsky space with some effective 
external field 

+ 72 o«,nMfl + Mn A**/J ~ ' + Mn/V./j) = °- ( 30 ) 
(2 Z 

3 Spontaneous breaking of symmetry for gravi- 
tons 



Let us consider vacuum solution <|30j) depending only on time. In quantum field 
theory this means dependence of vacuum on time. Then 



a 
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M/3,0,0 — M/3 l^l^nltf)- (31) 



a ^ a 4 



Taking into account constraints (|15p in variables /i} ~ n\ = Mi one 

gets the potential corresponding to (|3"Tj) as 



Write the field /i^ close to the minimum of the potential energy 

^=|^(0K(0), M £=/^(0)+$. (33) 

One must note that the condition (|33[) on (0) is the condition of minimal 
energy at same fixed moment t^. 

This is basic idia. Initied of dealing this time depended m, A we put the 
initial conditions at some to- Dusing these coordinations form the principle 
minimal energy at this moment. Surely m, A at this moment are numbers. 
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Take the solution (1551) as 



A4(0) = A*i(0) = 0, ^(0) = - M ^(0)=/xo = y^ r . (34) 

Then the Lagrangian 

1 a" n' 2 

can be written as 

l = ±04(0) + $)- n Mo) + &),n + f^(o) + ^)(/*£(o) + &)- 



a 
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^(^(oK(o) + $g + 2nZ(o)gy = 

+4/^^ + ^Ml? + (CI) 2 - 2£g)). 
Consider the quadratic in £2 terms 

i = ^ B 4» + ~ + ^ + fe2)2 " 2 ^ (35) 

,/2 „// 



Taking into account the equality A«o"~t = — anc ^ the gauge + £f = one 

a 4 a 



obtains the Lagrangian for gravitons 



L® = l® n £, n -^(&) 2 + &n (36) 

which is called by us the effective Lagrangian after the spontaneous breaking of 
symmetry. (46). Euler-Lagrange equations have the form 

On" On" 

^ + ^ = 0, ^ + ^1=0, g£ = 0, £5 = 0. (37) 



One sees that now in (|3T[) the sign of the mass squared is a correct one. The 
components £a>£f are massless while have the nonnegative mass squared 

2a" 

. The solutions for diagonal components £,,£2 (or = 1,2) can be 

a 

written as £ the Fourier integral 

= / dkic^Me 1 ^ + ctg k (v)e~ rkS ), (38) 
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And one has the equation 

g» + ^ + ^L. )gk = . (39) 
a 

This equation is free from the problem of the negative square of the effective 
2a" 

mass if > and one can construct the quantum theory of gravitons based 

a 

on the new vacuum state (f3"4"|) . One can notice that the vacuum expectation 
value of the field fj,(rf) is close to the scale factor. For the scale factor aijj) = rf 
one obtains 

Mo = a ^)\j^~~- 

One sees that for all p e [— 1;0) (if p S (0; 1) then to 2 = -2a" /a > 
and we leave vacuum as (1q = 0) the dynamical perturbation hp ^ 1 which 
contradicts the condition of the expansion of the curvature tensor into a series 
in this perturbation. The value p = — 1 corresponds to inflation, so we cannot 
deal the inflation model here while other situations can be considered. However 
this case must be considered separately and it is not studied in this paper. 

4 The Lagrange formalism for gravitons 

One can see from (|38M39p that gravitons in the expanding isotropic Universe are 
described by the effective scalar field $(x) with the Lagrangian 

L = V~g{${x), n d{xy n - ^Rd{x)d{x)), (40) 

where g = det(gtk) and R the scalar curvature. There is no term ^ because one 
deals with two polarizations. Euler-Lagrange equation for the field is £■[ = £| = 
£ = d ■ a. Look for solutions of (|39|) in the form (|38|) . The numbers C£,ci are 
changed on the operators c^,ci with commutation relations 

[c^c + ] = S(k-k'), [%,%] = [c+ c + ]=0. (41) 
The Fock vacuum state |in > is defined as 

c^|in >= 0, < in|in >= 1. 

Then for d(x) — —£,(x) one has 
a 

= ,„ ^ ^ / dk^Me 1 * 3 + c+g k ( V )e~ rkS ), (42) 



(27r) 3 a(?7) 
where gk(v) satisfy (|3"9"|) written as 



In" 

fl£ + w£(»7)fffc=0, w2fo) = — + k 2 (43) 

a 
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with initial conditions 



1 



9k(vo) = , = , g'k(Vo) = i\/uJ k (m)- (44) 

The condition for the wronskian 

9k(Vo)g'k* (Vo) ~ 9k(Vo)g'k(m) = ~2i (45) 



leads to existence of the full set of solutions of (j3"9")l in the sense of the indefinite 
scalar product 

(6,6) =ifdm fib 6)- (46) 
The Hamiltonian of the quantized field £(x) in the metric (|13[) has the form 

H(T,) = JdZ@ + ?+^£+t). (47) 

Putting the field d(x) from (gZJ) into (J32J) one obtains 



1 



#M = ~ / *?dkw k (ji)(E k (rj)(c+G k + c fe c+)- 



+ F fc ?+c++F fc *c fe c fc ), (48) 
where the coefficients -E^, are expressed through the solutions of the (|43[) 

= ^(\g' k \ 2 + M 2 ) Fk(v) = 7^(9' k 2 + gi)- (49) 

The corpuscular interpretation can be made in terms of creation and anni- 
hilation operators b k ,b k diagonalizing the Hamiltonian. If 

Ck = a* k (i])b k - Pk{v)bt 

then the Hamiltonian is 

^ I r°° ^ ^ 

H (^ = ^m k 2 dku; k (r 1 )(E k (r 1 )-l)(b+b k + b k b+). (50) 
n a (V) Jo 

The density of created particles and their energy density [4] can be found 
using formulas 

n (v) = ^rrr^ / k dk \^\ - ( 51 ) 

T Ol'j) JO 



10 



5 Some models of graviton creation 



Let us consider some matter filling the Universe with the equation of state p = je 
where p is pressure and e the energy density. One has for the homogeneous 
quasieuclidean isotropic Universe the equation [3] 

87tk 3a' 2 , . . 2 

—re = —7-, than am) = Cn , (52) 
c 4 a 4 

Let us take a(rf) — Crf and put it into (j4"3")l , then one obtains (However one 
must remember that due to our condition not all p can be used. One must have 



i 

n 

where 



"(,,) + (2p(p-l)- 3 + ) i 2 ).,(,,) = 0. (53) 



"' 2 = 2 *-^TJW' «»"-? 

So the results obtained for the scalar field in [lj are valid for gravitons for 
any scale factor with the scale factor of a given form. In [3] it was shown that for 
the density of created particles and the energy density defined by (f50"l - 15~1]) one 
gets convergent integrals. Let us calculate them. Putting the notation x = kf) 
one gets 

+ (1 + = 0. (54) 

Then the energy density of created particles due to (f5"4")l with is calculated 

as 

e{rf) =< 0|T °|0 >= 
x *dxu;(x)(-^(\^P-\ 2 + oj 2 (x)\g{x)\ 2 ) - 1) (55) 



7r 2 (a(?7)?7) 4 J 2lo(x) dx 

The solutions of (15411 are Bessel functions 



g(x) = Cl ^J{ l -V^m 2 ~,x) + C 2 ^Y{ l -^/l^^,x). 
Than 

2 , 1.5 ■ l(T 3 i? 3 / 2 , , . 

£M ~^W° Mm = a(r,)rf ' ° <m<4 ' (56) 

For the density of created particles in the unit volume one gets 



This integral is convergent [3]. For small m (0 < to < 0.5) 72(77) ~ R and for 
large to (to > 0.5) 72(77) ~ V^R- 
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Consider dust Universe with a(rj) = Crj 2 . Then 



4 • 10~ 3 , 

Sgrav. = ^ (c ) (58) 

For the background classical matter one has 

2- 10 84 cefc- 2 , _ 4N 

Smatt. = ^ ( C ) ( 59 ) 

So for Planckean time = 10~ 43 ce/c) the graviton energy density created 
from vacuum is some ten percent of the matter density while at the inflation 
time tin/ = 10 _36 cefc) it is only 10~ 14 of matter. These numbers are consistent 
with our approximation for the metric in perturbation theory. At the modern 
epoch one gets from (f58| that the energy flow from the time of the end of 
inflation U n f = 10 _36 cefc is 

Ssov.grav. = 0.5 ■ 10~ 12 (-^-j) (60) 

s • cm 

This can be compared with the flow from the Crab nebula [5] . One sees that 
it is much smaller. 

Skrab - 10" 8 ( *4) (61) 



6 Acknowledgements 

The authors are indebted to the participants of the A. A. Friedmann seminar 
of St. Petersburg for the discussions of the paper and to Ministry of Education 
and Science of Russia (grant RNP. 2. 1.6826) for financial support. 



References 

[1] Grishuk L. P. Relic gravitational waves and cosmology. UFN, vol. 175, 12, 
2005 (in Russian). 

[2] Lifshits E. M. Gravitational stabilities of the expanding world. JETP, 16, 
587 1946 (in Russian). 

[3] Landau L. D., Lifshits E. M. Field theory M.: Nauka, 1973. 

[4] Grib A. A., Mamayev S. C, Mostepanenko V. M. Vacuum quantum effects 
in strong fields. Friedmann Lab Publ., St. Petersburg, 1994. 

[5] Weinberg S. Gravitation and cosmology. M:Platon, 2000. 

[6] Dorofeev V. Yu. Expansion of Einstein equations of backgraund gravi- 
tational field, in Quantum theory and cosmology. Friedmann Lab Publ., 
St. Petersburg, 2009 (in Russian). 



12 



